WHAT IS CLAIMED IS: 



1. A method for determining signal parameters for one or more tones in an 
5 input signal, the method comprising: 

(a) receiving samples of the input signal, wherein the input signal includes the one 
or more tones; 

(b) operating on the samples to generate a transform array, wherein the transform 
array includes a positive frequency image and a negative frequency image for each of the 

10 one or more tones; 

(c) identifying frequency locations of one or more first magnitude peaks in the 
transform array; 

(d) computing a frequency estimate, ampUtude estimate and phase estimate for 
each of the one or more tones based on complex values of the transform array in a 

15 neighborhood of a corresponding one of the frequency locations; 

(e) correcting the complex values of the transform array in the frequency 
neighborhood of each frequency location based on the frequency estimates, amplitude 
estimates and phase estimates for the one or more tones; 

(f) computing an improved frequency estimate, improved amplitude estimate and 
20 improved phase estimate for each of the one or more tones based on the corrected 

complex values in the neighborhood of the corresponding frequency location; 

(g) storing the improved frequency estimates, improved amplitude estimates and 
improved phase estimates for the one or more tones. 

25 2. The method of claim 1, wherein (e) comprises: 

(el) estimating, for a first tone of said one or more tones, the complex 
contributions of positive and negative frequency images of tones other than said first tone 
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based on the corresponding frequency estimates, amplitude estimates and phase 
estimates; 

(e2) performing a pluraHty of summations resulting in a plurality of summation 
values, wherein each summation value includes one of the complex contributions for each 
5 of said positive and negative frequency images of tones other than said first tone; 

(e3) subtracting the sunmiation values from the corresponding complex values of 
the transform array in the frequency neighborhood of the frequency location 
corresponding to the first tone. 

10 3. The method of claim 2, wherein (e) further comprises: 

estimating, for the first tone, the complex contribution of the negative frequency 
image of the first tone based on the frequency estimate, ampUtude estimate and phase 
estimate of the first tone; 

wherein each sunmiation further includes one of the complex contributions of the 
15 negative frequency image of the first tone. 

4, The method of claim 2, wherein (e) further comprises performing (el), 
(e2) and (e3) repeatedly as the first tone runs through the one or more tones. 

20 5. The method of claim 1, 

wherein one or more of the positive frequency images and negative frequency 
images of the one or more tones at least partially overlap. 

6. The method of claim 1, wherein (e) operates to compensate the effects of 
25 the negative frequency images of the one or more tones on the complex values of the 
transform array in the frequency neighborhood of each frequency location. 
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7. The method of claim 1, wherein, for a first tone of said one or more tones, 
(e) operates to compensate the effects of the positive frequency images of tones other than 
said first tone on the complex values of the transform array in the frequency 
neighborhood of the frequency location corresponding to the first tone. 

5 

8. The method of claim 1, vv^herein said operating on the samples comprises: 
applying a window function to the samples to generate a windowed input signal; 

and 

computing a discrete Fourier transform of the windowed input signal to determine 
10 the transform array. 

9. The method of claim 1, wherein (d) comprises computing the frequency 
estimate and amplitude estimate for each of the one or more tones based on the 
magnitudes of the complex values of the transform array in the neighborhood of the 

1 5 corresponding frequency location. 

10. The method of claim 1, wherein (d) further comprises computing the phase 
estimate for each of the one or more tones based on the phase angle of one or more of the 
complex values of the transform array in the neighborhood of the corresponding 

20 frequency location. 

11. The method of claim 1, wherein (f) comprises computing the improved 
frequency estimate and the improved amplitude estimate for each of the one or more 
tones based on the magnitudes of the corrected complex values in the neighborhood of 

25 the corresponding frequency location. 

12. The method of claim 1, wherein (f) further comprises computing the 
improved phase estimate for each of the one or more tones based on the phase angle of 
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one or more of the corrected complex values in the neighborhood of the corresponding 
frequency location. 

13. The method of claim 1, further comprising: 
repeating (e) and (f) one or more times before performing (g); 

wherein the improved frequency estimates, improved amplitude estimates and 
improved phase estimates from a current iteration of (f) are used as the frequency 
estimates, amplitude estimates and phase estimates in a next iteration of (e). 

14. The method of claim 1, wherein, for a first tone of said one or more tones, 
(e) comprises: 

(el) computing a first shifted window value W(k,j) at a bin index value k for each 
negative frequency image Nj based on the frequency estimate for the corresponding tone 

(e2) multiplying each first shifted window value W(k,j) by the amplitude estimate 
for corresponding tone xj to generate a first product G(k,j); 

(e3) phase shifting each first product G(kj) by an angle which depends on the 
phase estimate of the corresponding tone xj; 

(e4) computing a summation value D(k), wherein the sunamation value D(k) 
includes the first product G(k,j) for each tone xj of the one or more tones; 

(e5) subtracting the sunmiation value D(k) from the corresponding complex value 
of the transform array in the frequency neighborhood of the frequency location 
corresponding to first tone. 

15. The method of claim 14, wherein, for the first tone of said one or more 
tones, (e) further comprises: 
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computing a second shifted window value T(kj) at the bin index value k, for each 
positive frequency image Pj except the positive frequency image P] corresponding to the 
first tone, based on the frequency estimate for the corresponding tone Xj; 

multiplying each second shifted window value T(k,j) by the amplitude estimate 
for corresponding tone Xj to generate a second product G(k,j); 

phase shifting each second product G(k,j) by an angle which depends on the phase 
estimate of the corresponding tone Xj; 

wherein the sunomation value D(k) includes the second product G(k,j) for each 
tone Xj of the one or more tones other than the first tone. 

16. The method of claim 14, further comprising: 

repeating (el) through (e5) for each bin index value k in the frequency 
neighborhood of the frequency location corresponding to the first tone. 

17. The method of claim 1, wherein (c) comprises: 

determining magnitude peaks in the magnitude spectrum of the transform array 
which exceed a magnitude threshold in a range of positive frequencies; and 

computing a center frequency for each of the threshold-exceeding magnitude 

peaks. 

18. The method of claim 1, further comprising: 

transmitting an indication of the improved frequency estimates, improved 
ampUtude estimates and improved phase estimates for the one or more tones to an output 
device. 

19. A system for determining signal parameters for one or more tones in an 
input signal, the system comprising: 



Atty, Dkt. No,: 5150-42800 



Page 36 



Conley, Rose & Tayon, P.C. 



an input for receiving samples of the input signal, wherein the input signal 
includes the one or more tones; 
an output device; 

a memory configured to store program instructions; 
5 a processor configured to read the program instructions from the memory and to 

execute the program instructions, wherein, in response to execution of the program 
instructions, the processor is operable to; 

(a) receive samples of the input signal, wherein the input signal includes the one 
or more tones; 

10 (b) operate on the samples to generate a transform array, wherein the transform 

array includes a positive frequency image and a negative frequency image for each of the 
one or more tones; 

(c) identify frequency locations of one or more first magnitude peaks in the 
transform array; 

15 (d) compute a frequency estimate, amplitude estimate and phase estimate for each 

of the one or more tones based on complex values of the transform array in a 
neighborhood of a corresponding one of the frequency locations; 

(e) correct the complex values of the transform array in the frequency 
neighborhood of each frequency location based on the frequency estimates, amplitude 

20 estimates and phase estimates for the one or more tones; and 

(f) compute an improved frequency estimate, improved amplitude estimate and 
improved phase estimate for each of the one or more tones based on the corrected 
complex values in the neighborhood of the corresponding frequency location; and 

(g) transmit an indication of the improved frequency estimates, improved 
25 amplitude estimates and improved phase estimates for the one or more tones to an output 

device. 

20. The system of claim 19, 
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wherein one or more of the positive frequency images and negative frequency 
images of the one or more tones at least partially overlap. 

21. The system of claim 19, wherein (e) operates to compensate the effects of 
5 the negative frequency images of the one or more tones on the complex values of the 

transform array in the frequency neighborhood of each frequency location. 

22. The system of claim 19, wherein, for a first tone of said one or more tones, 
(e) operates to compensate the effects of the positive frequency images of tones other than 

10 said first tone on the complex values of the transform array in the frequency 
neighborhood of the frequency location corresponding to the first tone. 

23. The system of claim 19, wherein, in operating on the samples, the 
processor is operable to: 

15 ^ppiy a window function to the samples to generate a windowed input signal; and 

compute a discrete Fourier transform of the windowed input signal to determine 
the transform array. 

24. The system of claim 19, wherein, in response to the execution of the 
20 program instruction, the processor is operable to compute the frequency estimate and 

amplitude estimate for each of the one or more tones based on the magnitudes of the 
complex values of the transform array in the neighborhood of the corresponding 
frequency location. 

25 25. The system of claim 19, wherein, in response to the execution of the 

program instructions, the processor is operable to compute the phase estimate for each of 
the one or more tones based on the phase angle of one or more of the complex values of 
the transform array in the neighborhood of the corresponding frequency location. 
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26. The system of claim 19, wherein, in response to the execution of the 
program instructions, the processor is operable to compute the improved frequency 
estimate and the improved amplitude estimate for each of the one or more tones based on 

5 the magnitudes of the corrected complex values in the neighborhood of the corresponding 
frequency location, 

27. The system of claim 19, wherein in response to the execution of the 
program instructions, the processor is operable to compute the improved phase estimate 

10 for each of the one or more tones based on the phase angle of one or more of the 
corrected complex values in the neighborhood of the corresponding frequency location. 

28. The system of claim 19, wherein the output device is a display device. 

15 

29. A method for determining signal parameters for one or more tones in an 
input signal, the method comprising: 

(a) receiving samples of the input signal, wherein the input signal includes the one 
or more tones; 

20 (b) operating on the samples to generate a transform array, wherein the transform 

array includes a positive frequency image and a negative frequency image for each of the 
one or more tones; 

(c) identifying frequency locations of one or more first magnitude peaks in the 
transform array; 

25 (d) computing one or more of a frequency estimate, amplitude estimate and phase 

estimate for each of the one or more tones based on complex values of the transform array 
in a neighborhood of a corresponding one of the frequency locations; 
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(e) correcting the complex values of the transform array in the frequency 
neighborhood of each frequency location based on one or more of the frequency 
estimates, amplitude estimates and phase estimates for the one or more tones; 

(f) computing one or more of an improved frequency estimate, improved 
5 amplitude estimate and improved phase estimate for each of the one or more tones based 

on the corrected complex values in the neighborhood of the corresponding frequency 
location; 

(g) storing the one or more of the improved frequency estimates, improved 
amplitude estimates and improved phase estimates for the one or more tones. 

10 

30. The method of claim 29, further comprising: 

transmitting an indication of one or more of the improved frequency estimates, 
improved amplitude estimates and improved phase estimates for the one or more tones to 
an output device. 

15 

31. A memory medium comprising program instructions for determining 
signal parameters for one or more tones in an input signal, vv^herein the program 
instructions are executable by one or more processors to implement: 

20 (a) receiving samples of the input signal, wherein the input signal includes the one 

or more tones; 

(b) operating on the samples to generate a transform array, v^herein the transform 
array includes a positive frequency image and a negative frequency image for each of the 
one or more tones; 

25 (c) identifying frequency locations of one or more first magnitude peaks in the 

transform array; 
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(d) computing a frequency estimate, amplitude estimate and phase estimate for 
each of the one or more tones based on complex values of the transform array in a 
neighborhood of a corresponding one of the frequency locations; 

(e) correcting the complex values of the transform array in the frequency 
5 neighborhood of each frequency location based on the frequency estimates, amphtude 

estimates and phase estimates for the one or more tones; 

(f) computing an improved frequency estimate, improved amphtude estimate and 
improved phase estimate for each of the one or more tones based on the corrected 
complex values in the neighborhood of the corresponding frequency location; 

10 (g) transmitting an indication of the improved frequency estimates, improved 

amplitude estimates and improved phase estimates for the one or more tones to an output 
device. 

32. The memory medium of claim 31, wherein (e) operates to compensate the 
15 effects of the negative frequency images of the one or more tones on the complex values 

of the transform array in the frequency neighborhood of each frequency location. 

33. The memory medium of claim 31, wherein, for a first tone of said one or 
more tones, (e) operates to compensate the effects of the positive frequency images of 

20 tones other than said first tone on the complex values of the transform array in the 
frequency neighborhood of the frequency location corresponding to the first tone. 

34. The memory medium of claim 31, wherein (d) further comprises 
computing the phase estimate for each of the one or more tones based on the phase angle 

25 of one or more of the complex values of the transform array in the neighborhood of the 
corresponding frequency location. 
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35. The memory medium of claim 31, wherein (f) further comprises 
computing the improved phase estimate for each of the one or more tones based on the 
phase angle of one or more of the corrected complex values in the neighborhood of the 
corresponding frequency location. 

5 

36. The memory medium of claim 3 1 , further comprising: 
repeating (e) and (f) one or more times before performing (g); 

wherein the improved frequency estimates, improved amplitude estimates and 
improved phase estimates from a current iteration of (f) are used as the frequency 
10 estimates, amplitude estimates and phase estimates in a next iteration of (e). 

37. A method for determining signal parameters for a plurality of tones in an 
input signal, the method comprising: 

(a) receiving samples of the input signal, wherein the input signal includes the 
15 plurality of tones; 

(b) operating on the samples to generate a transform array comprising complex 

values; 

(c) computing a frequency estimate and an amplitude estimate for each of the 
plurality of tones based a corresponding first magnitude peak in the magnitude spectrum 

20 of the transform array; 

(d) computing a phase estimate for each of the pluraUty of tones based on the 
phase of at least one of the complex values of the transform array in a frequency 
neighborhood of the corresponding first magnitude peak; 

(e) correcting the complex values of the transform array in the frequency 
25 neighborhood of each first magnitude peak based on the corresponding frequency 

estimates, amplitude estimates and phase estimates of the plurality of tones; 
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(f) computing an improved frequency estimate and improved amplitude estimate 
for each of the plurality of tones based on a corresponding second magnitude peak of the 
corrected complex values; 

(g) computing an improved phase estimate for each of the plurality of tones based 
5 on the phase of at least one of the corrected complex values in the frequency 

neighborhood of the corresponding second magnitude peak; 

(h) transmitting as output the improved frequency estimate, improved amplitude 
estimate and improved phase estimate. 

10 38. The method of claim 37 further comprising: 

repeating (e) through (g) until a convergence condition based on one or more of 
improved frequency estimate, improved amplitude estimate, and improved phase 
estimate. 

15 39. The method of claim 38 wherein the improved frequency estimates, the 

improved amplitude estimates, and the improved phase estimates determined in a current 
iteration of (f) and (g) as used as the frequency estimates, amplitude estimates and phase 
estimates in a next iteration of (e). 

20 40. The method of claim 37, 

wherein the transform array comprises an additive combination of a positive 
frequency image and a negative frequency image corresponding to each of the pluraUty of 
tones; 

wherein (e) comprises: 

25 estimating, for a first tone, the complex contributions of positive and 

negative frequency images of tones other than said first tone based on the corresponding 
frequency estimates, amplitude estimates and phase estimates; 
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performing a plurality of summations resulting in a plurality of summation 
values, wherein each summation value includes one of the complex contributions for each 
of said positive and negative frequency images of tones other than said first tone; 

subtracting the summation values from the corresponding complex values 
5 of the transform array in the frequency neighborhood of the first magnitude peak 
corresponding to the first tone. 

41. The method of claim 40, wherein (e) further comprises: 
estimating, for the first tone, the complex contribution of the negative frequency 
10 image of the first tone based on the frequency estimate, amplitude estimate and phase 

estimate of the first tone; 

wherein each sunmiation further includes one of the complex contributions of the 

negative frequency image of the first tone. 

15 42. The method of claim 37, wherein said operating on the samples comprises: 

applying a window function to the samples to generate a windowed input signal; 

and 

computing a discrete Fourier transform of the windowed input signal to determine 
the transform array. 

20 

43. A method for determining signal parameters for a tone comprised within 
an input signal, the method comprising: 

(a) receiving samples of the input signal; 

(b) operating on the samples to generate a transform array comprising complex 

25 values; 

(c) computing a frequency estimate and an amplitude estimate for the tone based a 
first magnitude peak in the magnitude spectrum of the transform array; 
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(d) computing a phase estimate for the tone based on the phase of at least one of 
the complex values of the transform array in a frequency neighborhood of the first 
magnitude peak; 

(e) correcting the complex values of the transform array in the frequency 
5 neighborhood of the first magnitude peak based on the frequency estimate, amplitude 

estimate and phase estimate; 

(f) computing an improved frequency estimate and improved amplitude estimate 
for the tone based on a second magnitude peak of the corrected complex values; 

(g) computing an improved phase estimate for the tone based on the phase of at 
10 least one of the corrected complex values in the frequency neighborhood of the 

corresponding second magnitude peak; 

(h) transmitting as output the improved frequency estimate, improved amplitude 
estimate and improved phase estimate. 

15 44. The method of claim 43 further comprising: 

repeating (e) through (g) until a convergence condition based on one or more of 
improved frequency estimate, improved amplitude estimate, and improved phase 
estimate. 

20 



25 
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